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ABSTRACT 

Starting from an unitary effective Lagrangian model for the meson-nucleon scattering 
developed in we come to a unified description of both meson scattering and photon- 
induced reactions on the nucleon. To this end the photon is added perturbatively, yielding 
both Compton scattering and meson photoproduction amplitudes. In a simultaneous fit to 
all available data the parameters of the nucleon resonances are extracted. We find that a 
global fit to the data of the various channels involving the final states 7A^, vriV, irnN, rjN 
and KA is possible. Especially in the eta photoproduction a readjustment of the masses and 
widths found in the fits to hadronic reactions alone |l|] is necessary to describe the data. Only 
for the £'13(1520) we find a possible disagreement for the helicity couplings extracted using 
the combined dataset and pion photoproduction multipoles alone. The model-dependence 
introduced by the restoration of gauge invariance is discussed and found to be significant 
mainly for resonances with small helicity couplings. 
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I. INTRODUCTION 



Photon- and meson-induced reactions on the nucleon are the main source of information 
about the nucleon resonance spectrum. From the knowledge of the possible excitations of 
the nucleon and their properties one hopes to extract information about the structure of 
the nucleon. To this end one needs models that allow to determine the masses and partial 
decay widths of the resonances. Because of the constraint of unitarity the different reaction 
channels can in principle not be treated separately, but have to be taken into account 
simultaneously. 

For the purely hadronic reactions models based on unitarity and analyticity have been 
widely used 0-0. Using an ansatz developed by the Carnegie- Mellon Berkeley or CMB 
group, Dytman et al. for example have recently extracted resonance parameters from a fit 
to the TxN 7iN, niiN, f]N data [||. 

On the other hand, in calculations of meson photoproduction, effective Lagrangian mod- 
els are the main tool for investigations p|-p!2[]. In these models the important constraint of 
gauge invariance can be easily implemented on the operator level. However, unitarity has 
only been fulfilled in a few calculations of pion photoproduction in the A-region P JT0|JT3|JT^ 



and eta photoproduction [|T^]. Lately, also a description of Compton scattering in this 
framework has been put forward [|r^. 

With the availability of high-precision data from various accelerator facilities like MAMI, 
ELSA, GRAAL and TJNAF, there is also an urgent need to improve the models for meson 
photoproduction. The most important ingredient for these improved models is the dynamical 
treatment of the meson rescattering in the same framework as the initial photoproduction 
reaction. To this end a description of the purely hadronic reactions within the effective 
Lagrangian approach is necessary. 

As a first step, we formulated a model for these hadronic reactions employing the K- 
matrix approximation The interaction potential is described solely in terms of s-, u- and 
t-channel Feynman diagrams. We have shown that a reliable extraction of the resonance 
parameters is possible in this model and that the effective Lagrangian approach reduces the 
number of free parameters considerably, since the non-resonant background is generated by 
a few couplings only. 

Analyticity is not guaranteed in this approach, but there are estimates about the quality 
of the i^'-matrix calculation as compared to other approximations . These indicate that 
the final resonance parameters extracted in different approximations are very similar. 

The aim of this paper is now, to extend this model to include also photon-induced 
reactions. This would allow to benefit from the very accurate data for such reactions that 
already are or will be available in the near future. To this end we will at first shortly discuss 
the model used and the results of the fits to hadronic data. The extension to photon-induced 
reactions and the database available for the various reactions will then be addressed. After 
that, the results of the fits to the combined data are presented and discussed. 

An important new feature of this analysis is the extraction of electromagnetic coupling 
constants from a combined fit to meson photoproduction and Compton scattering data. 
Apart from the t-channel contributions, the latter is determined only by the electromag- 
netic couplings and not by a product of strong and electromagnetic couplings as in meson 
photoproduction. Compton scattering data should thus provide an important constraint on 
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the extraction of electromagnetic coupling constants. We will also use a comparison with a 
dispersion theoretical analysis of these data to obtain information about the quality of the 
i^-matrix approximation. 



II. DESCRIPTION OF THE MODEL AND ITS APPLICATION TO THE 
PURELY HADRONIC REACTIONS 

For easier reference we review in this section briefly the treatment of the hadronic chan- 
nels in Ijl]]. This then forms the basis of our inclusion of photon- induced reaction channels. 

A. The JC-matrix approximation 

The Bethe-Salpeter equations encountered in meson-nucleon scattering 

T = V + VGbsT (1) 
can always be decomposed into two coupled equations p|,p!6|: 

K = V + VRe{GBs)K 

T = K~ iKlm{GBs)T. (2) 

Here Gbs is the Bethe-Salpeter propagator that describes the intermediate propagation of 
the meson and the nucleon. 

In (ID the second equation can easily be solved, since the imaginary part of Gbs is always 
proportional to 5-functions that place the meson and nucleon on its mass-shell. The real part 
of Gbs in the first equation in (H) amounts to a principle- value integral and is, therefore, 
much harder to treat. 

The ii'-matrix approximation now consists of neglecting this principle-value integral in 
the first equation and thus using = in the second one. This leads to 



[Tk] = [V + iVTK] 



V 



iV 



(3) 



where the brackets [■ ■ ■] indicate that V and are n x n-matrices [n being the number of 
asymptotic channels taken into account) and that (0) is a matrix equation. 

Obviously, Tk as calculated from (|^) does not necessarily fulfill a dispersion relation, 
as does the full T from (^. Therefore, the iiT- matrix approximation does not guarantee 
analyticity by construction. However, Pearce and Jennings 0] have shown that in nN- 
scattering the contributions from Iie{GBs) are of minor importance, since the corresponding 
principle- value integral is reduced by a very soft cutoff needed to regularize (^. Furthermore, 
Surya and Gross estimated that in this process the error made by putting the pion onshell 
is of the order of {Et^ — m^)£'^/m^ and can therefore be neglected for small pion energies. 
Based on these studies, it seems that the main contributions in the Bethe-Salpeter equation 
(0) come from the rescattering with both intermediate particles onshell. This part is correctly 
taken into account in the iiT-matrix approximation. 
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B. Contributions to the potential V 



As asymptotic states vrA^, miN, rjN and KA have been taken into account. A detailed 
description of the database used in the fits is given in Neglected so far are final states 
such as KT, and uN. This has been done because either the coupling of resonances to this 
channel is known to be small {KT,, [0]), or because only one resonance is known that might 
have a significant decay into this channel {ujN, p). 

Also the pA^-channel is not contained in our analysis. To make at least partly up for this 
deficiency, we parameterize the vrvrA^-decay by the coupling to a scalar, isovector (^-meson 
P| JT^JTB|JTE| ] with mass = 2m.,^. This ensures that the important vrvrA^-decay is taken into 
account, but at the same time the model is kept as simple as possible. A description of the 
TTirN final state in terms of two particle intermediate states like ttA, 7rPii(1440) and pN is 
currently under way. 

The potential V is now calculated from the interaction Lagrangians collected in App. 
Taken into account are contributions from the nucleon Born terms, t-channel exchanges of 
p, ao and K* and resonance contributions in the s- and u-channels. 

In this work we limit ourselves to partial waves with spin ^ and |. This is because only 
for these the Lagrangians can be given in an unambiguous way Furthermore, we 



restrict the energy range to ^/s < 1.9 GeV because for higher energies additional decay 
channels might be important. 



C. Form factors 

In order to investigate the dependence of the resonance parameters on the specific choice 
of the hadronic form factors, in [Q] fits using combinations of the three basic forms 

A'' 



A4 + (g2 _ ^2)2 



Fe(.^)^exp(-i^!^) 

^ A4+(g2_(g2^^2/2))2 ^ ^ 

have been carried out (m denotes the mass of the propagating particle, q its four-momentum 
and qf is the value of at the kinematical threshold in the t-channel). The following 
restrictions were imposed to limit the number of free cutoff parameters: 

• the same functional form F and cutoff A^r is used in all vertices ttNN, rjNN and 
KNA, 

• for all resonances we take the same F as for the nucleon, but different values A1/2 and 
A3/2 for the cutoffs for spin-| and spin-| resonances, 

• in all t-channel diagrams the same F and At are used. 
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D. Results for the hadronic reaction channels 



In this section we briefly recapitulate the main results of the fits to the hadronic reactions 
TiN — i> ttA^, ttN —y TTTiN, irN — >■ riN and nN — >■ KA. The reader is referred to [|[] for a more 
detailed discussion. 

• A qualitative and quantitative description of all hadronic data is achieved. 

• The resonance parameters we find are in good agreement with the values obtained by 
other groups [g|,|,0,|T8| and with the PDG-values ^ 



• The IT p ^ r]n data are not good enough to determine the r^A^-branching fractions of 
the nucleon resonances very accurately. Nevertheless, we find non-vanishing couplings 
for a few resonances (^ii(1535), Pii(1710), Pi3(1720) and L)i3(1520)). 

• Only two resonances (5*11(1650) and Pii(1710)) exhibit sizeable couplings to the A'A- 
channel. For higher energies the tt~p A^^A-reaction is dominated by the t-channel 
A'*-exchange. 

• There are deviations from the vrA^ data for the P/3- and D/3-partial waves. Below the 
resonance positions we underestimate the amplitudes, whereas for energies above the 
resonance the fits give too large amplitudes. This indicates that a fit with form factors 
that are asymmetric around the resonance position rriR (instead of our choices Pp^g) 
might lead to a better description of the data. 

• Already for energies above 1.6 GeV we find that the t-channel contributions are not 
adequately described by the corresponding Feynman diagrams anymore. It seems 
that stronger modifications than those from any of the form factors Pp,e,i are needed. 
A smooth transition to a Regge-like behavior in this energy range would probably 
improve the quality of the fits for the highest energies. 

Because of this, in the fits the ttNN- and pA^A^-coupling constants are driven below 
their usual values. This in turn also leads to a smaller P33(1232) mass (1.229 GeV) 
and width (110 - 113 MeV). 

A comparison with the resonance parameters obtained from other unitary and analytic 



analyses, as for example those listed in [17|, shows that the hadronic reactions can be 



described rather well in this effective Lagrangian approach employing the AT-matrix approx- 
imation. 



III. EXTENSION TO PHOTON-INDUCED REACTIONS 

With the inclusion of the 7A^ final state a model can now be constructed that combines 
the advantages of describing the electromagnetic interactions using effective Lagrangians 
with a dynamical treatment of the rescattering. In principle this extension is straightforward, 
it mainly consists of enlarging the matrices [V] and [T] to take into account the new final 
state. Eq. @j then gives the unitarized amplitudes for meson photoproduction and Compton 
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scattering. However, there are four important points that make the combined treatment of 
all possible reaction channels technically more involved: 

1. The photon can induce electric and magnetic transitions. Therefore, in the case of 
spin- 1 resonances two amplitudes (called multipoles in photoproduction, instead of 



partial waves as in the hadronic reactions) have to be taken into account |]T6 

2. Furthermore, the interaction of the photon with the nucleon can be split up into an 
isoscalar and an isovector part. For e.g. pion photoproduction this leads to three 
different isospin amplitudes (T°, T^/^ and T^/^), instead of two (T^/^ and T^/'^) as in 
vrA^-scattering |1T6| . 

3. In the multipole-decompositions of photoproduction data any influence of Compton 
scattering is usually neglected [plj] . 

4. The Compton amplitudes cannot be fully isospin-decomposed using experimental data 
alone. Since only two physical processes (7^ and 7n) exist, only two of the four 
amplitudes can be extracted. Therefore, the rescattering contributions are normally 
calculated in a basis using physical states (e.g. vr~p, vr^p, n^n, vr^n) and not pure 
isospin states ||21[ . 



The first point can easily be taken into account by introducing two new final states 
{'jN)e and {'yN)M, where the index denotes the type of electromagnetic transition. The 
isoscalar /isovector nature of the photon can be treated in a similar way. The usual isospin 
decomposition of pion photoproduction, for example, is given by |1^,|T6|: 



(7r,|T.,|7) = T,T^, + \r,r,Tl/' + (5,3 - \r,r,)Tl/'. (5) 

Here T^^ contains the amplitude for isoscalar photons, whereas T^/^'^/^ are the amplitudes 
for total isospin I = \i\ induced by isovector photons. The important point to note is, that 
rescattering of and T^/^ takes place via the I = \ part of the hadronic reaction channels 
(e.g. T^l"^), whereas for T^/^ only T^l^ contributes. The / = | amplitudes can, therefore, 
be treated in the usual way. For the / = | sector a further splitting of the final states into 

(7A^)^, (7A^)l/^ {lN)% and (7A^)if needs to be introduced. 

The third and the fourth point amount to neglecting photon-rescattering contributions 
when extracting the photoproduction amplitudes from the data. Only a direct term ~ V^^ 
is then present in Eq. (^. In our analysis of those amplitudes this can only be taken into 
account by using the form T = V + iTV and not T = V/{1 — iV), since in the latter case 
direct and rescattering terms cannot be disentangled. For the photoproduction of a scalar 
meson (p this schematically leads to 

^<^7 = + i X! ^w'K/j'V A ^ 0' 2 ' 2 ' ^ 2 ' 2 ' 
For Compton scattering we end up with: 

T,T = ^7T + iE^7cK7- (7) 
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Here the sum runs over all physical intermediate states c (e.g. vr^p, vr+n, ... for 7p — >■ 
7p). Neglected are in both cases contributions from the photon rescattering, since they 
are suppressed by an additional factor e^. Because of the same suppression we also do not 
take electromagnetic corrections to the hadronic channel into account. We have checked this 
approximation by also calculating the photon-rescattering contributions in the A-region and 
found them to be negligible. 



A. Background contributions and resonance couplings 

The contributions to the potential V in the case of photon-induced reactions come from 
Bremsstrahlung of asymptotic particles {N, ir, K) and electromagnetic decays (e.g. nucleon 
resonances and vector mesons). The Bremsstrahlung leads to Born diagrams in the s-, u- 
and t-channels; decays of nucleon resonances give contributions in the s- and u-channels, 
whereas meson decays enter as t-channel diagrams. 

To calculate the different contributions, we first of all need to specify the couplings to 
the photon. Since the Lagrangians have already been discussed many times [10,^,11,12,22 
we limit ourselves to a short summary. 

For the nucleon (N, A and S) and the scalar mesons (denoted by if) we have: 



CyNN 



-eN 



eNnA- 



N 



a 



' 4:171 N 

-e[cpxid,cp)].^A'^ 



2m 



N 



(8) 

(9) 
(10) 

(11) 



The magnetic transitions moments are given by = —0.613 and k^, = 1.610 0,^. In the 
case of pseudovector (^A^A^-coupling, C^^^nn gives rise to the so-called seagull- or 4-point- 
diagram. 

For the decay of scalar {(p) and vector {v) mesons we use [F^'^ = d'^A^ — d^A'^ , v^'^ = 



C 



77V 



e'^e,„,.,^F^%d'^)A'' 



2m,, 



(12) 



with the couplings extracted from the corresponding decay widths |[T7[| : 



-0.044, g^^yrj = 0.167 

0.131, c/^^±p± = 0.103, 

0.313, c/^^^ = 0.329 

0.631, c/7^±i^.,± = 0.415. 



1.020 



For the spin-| resonances only a magnetic transition is possible: 



(13) 
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C^NR,^, = eRgi^NF^" + h.c, (14) 



where R is the resonance spinor. The operator T^^, is given by 750"^,^ for odd-parity resonances 
and by cr^j, otherwise. In the spin-| case two couphngs can contribute: 

c^NRy, = -^i?"e„^(zi)7.riVF-'^ + 

-^R-e^^(z2)Tid,N)F^f^ + h.c, (15) 

4:171 

Qa,,{z) = g^^~^{l + 2z)7„7^. 

Here, the operator F is either 1 or 75 for odd and even parity resonances, respectively. 

The initial values for the electromagnetic couplings of the resonances are calculated from 
the helicity couplings given by the Particle-Data-Group []T^. The connections between the 
gi^2 used here and the Ai/2,3/2 are as follows p^j2l 

Spin- 1 

3 




Spin-I : Ay, = -i^^l -^^^^^ {±g, + ^{run T m^)92) 

{91- -r—{mRTmN)g2)- (16) 

Here S,r denotes the phase at the nNR-veTtex. These relations can easily be inverted to 
extract the gi^2- 

In all these couplings we have to account for the isoscalar/isovector nature of the photon. 
This is done in the case of J = | resonances by using gi^2 = ^712 + '^39i2 (P^ and (|T5|). 
Note that / = | resonances only couple to isovector photons and therefore have the same 
coupling 9i^2T3 to •yp and •jti. 

The formulae needed for the multipole decomposition of photoproduction and Compton 
scattering are somewhat lengthy and have been collected in App. p. There we also list the 
connections of the amplitudes to the various observables. 



B. Form factors and gauge invariance 

As we have already stated in the introduction, an effective Lagrangian model allows to 
address the question of gauge invariance on a fundamental level. Unfortunately, there is no 
unique way to restore gauge invariance once hadronic form factors have been introduced. 
This ambiguity affects only the Born terms, since all transition vertices fulfill gauge in- 
variance by construction [|]. That is, having a vertex function F^^^ constructed from the 
corresponding Lagrangian, we always have 



8 



where is the photon-momentum. The main reason for this is that a R ^ N'y transi- 
tion cannot be derived from the hadronic Lagrangians by minimal substitution, but has to 
be constructed by hand. Therefore, all corresponding vertex functions have to fulfill (|1^ 
separately, in order to preserve gauge invariance. Because of this, we have also introduced 
independent form factors at the electromagnetic vertices of the nucleon resonances. These 
form factors are taken to have the same analytic form as the hadronic ones (cf. Eq. @), 
only the cutoffs A are chosen independently. 

For vertices derived through minimal substitution, only the sum of all contributing (s-, 
u- and t-channel and 4-point) diagrams needs to fulfill a constraint similar to (p!7|): 



E 

=s,u,tA 



0. 



(18) 



Here Aii^^^i denotes the contribution of the ?'th Born diagram to the photoproduction am- 
plitude. 



In the case of pion photoproduction the vertex functions are F^^^, 



r^^^ and V^^NN and 

lead to four possible diagrams: s-, u- and t-channel and 4-point. Now the momenta at the 
hadronic vertices are different in all diagrams. Since, correspondingly, the form factors have 
different values for each diagram, all terms in ([1^) are weighted differently once hadronic 
form factors are introduced. Hence, the inclusion of any of the factors -Fp e,t will results in a 
violation of gauge invariance. 

The first model that tries to solve this problem was proposed by Ohta 



In his ansatz 

it is assumed that the form factor F can be separately Taylor-expanded with respect to 
the momenta p p' and q. After minimal substitution in all three momenta, the resulting 
expressions can then be resummed in closed form. 

To illustrate this in some detail, we focus on the s-channel Born diagram. For our choice 
of hadronic form factors the amplitude A^^^ is given by (suppressing factors of i): 



u{p')gT,NNF{s)'-f^ 



, (jf^ + + niN 
{p + - 171% 



10" 



2mM 



u{p). 



(19) 



Using now Ohta's prescription we obtain an additional counter term M. 
rearrangements the sum of both can be expressed as: 



Ohta 



After some 



Ohta 



+ 



u{p)g-KNNibl 



Ohta 

{p + ky - m%' 



ej^ulp) + Terms 



(20) 



From this it is clear that the influence of the form factor F on the coupling to the charge 
of the nucleon has been fully removed by the Ohta-prescription. Only the coupling to the 
magnetic moment (which is gauge invariant by itself) is affected by the introduction of F 
at the hadronic vertex. The same also holds for the u- and t-channel diagrams. 



Haberzettl has argued against this procedure because four-momentum conservation 



connects p, p' and q and the form factor F is needed at an unphysical point F(m 
To incorporate both points, he has developed an alternative method that leads to different 
expressions for the final amplitude. 
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If we clioose, for example, p and p' as the independent variables, we can remove the 
g-dependence by using q = p' — p. After this, we get one counter-term less than in Ohta's 
method. The net result for the electric s-channel contribution is in this case given by: 



Haberzettl ■^■y^,s 



■ F{t). (21) 

Ohta ^ ^ ^ ^ 



Instead of removing the influence of F altogether, we now have replaced the form factor 
used in the bare Born diagram with the one from the t-channel exchange. Since the same 
happens for the other diagrams as well, we have an overall factor F{t) multiplying the charge 
contributions from all four Born diagrams. Of course, the choice of independent variables 
is arbitrary, so that in its most general form the amplitude (|21|) contains, instead of F{t), 
a common form factor F = F{s, u, t). This degree of freedom has for example been used in 
the calculation of Nozawa et al. 0, where a form factor F = F{s) was taken into account. 
Following Haberzettl [^, we use here an ansatz of the form 

F{s, u, t) = aiF(s) + a2F{u) + asFit), (22) 

with the additonal constraint ai + 02 + as = 1; in [^] the coefficients were fitted in kaon 
photoproduction. 

To investigate the dependence of the electromagnetic couplings of the resonances on the 
gauge procedure, we have performed fits using both Ohta's method and the prescription of 
Haberzettl. Since this is mainly an exploratory study we have in the latter case adopted 
the 'democratic' choice and taken the a^'s to be equal. For pion photoproduction this leads 
to the ansatz F{s,u,t) = l/3(F(s) + F{u) + F(t)) (p2|). In the case of eta and kaon 
photoproduction, however, we use a slightly different form factor F because in both cases 
the Born terms consist of a smaller number of diagrams. For the eta case we don't have a 
t-channel contribution, whereas for the kaon photoproduction there is no w-channel diagram 
involving a coupling to the charge of the particles. Therefore, in our calculation the ansatz 

F^{s,u,t) = |F(s) + |F(m) + iF(t) 

F,{s,u) = |F(s) + iF(M) (23) 
FK{s,t) = |F(s) + |F(t), 

was used for the different form factors. This prescription has been chosen to avoid additional 
free parameters. But we want to stress again that (^) and (^) can only be motivated by 
the structure of the hadronic vertices present in the photoproduction diagrams. Only in 
a microscopic model for the electromagnetic vertex the form factor could be determined 
unambiguously. 

Finally, one could think about introducing form factors at the electromagnetic vertices 
of the nucleon and the pion and kaon as well. For the 'jNN case this would lead to the 
following vertex: 



Knn = u{p')e 



F,{p\p'\k')^^ + F,{p\p'\e) 



u{p) (24) 



For real photons gauge invariance demands Fi(p^,p'^,0) = 1. In addition, for F2 we have 
0) = kn- Besides this, we have no further constraint on F2 for the case of 
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meson photoproduction. In the case of Compton scattering it can be shown |2^, however, 
that additional contributions are needed to restore gauge invariance, once a form factor 
F2 has been introduced. Unfortunately, also these contributions cannot be determined in 
an unambiguous way. To avoid this additional uncertainty, we have chosen not to include 
electromagnetic form factors for the nucleon and scalar mesons in our calculation. 



C. Calculation of photon-induced reactions 

The calculation of photon-induced reactions is now carried out as follows: i) calculate 
the potential \y] from all contributing Feynman diagrams; ii) invert the hadronic submatrix 
to give the T-matrix {Thadr] = \^hadr/{)- — iV^adr)]; iii) unitarize the meson photoproduction 
with the help of (^); iv) finally, calculate the T-matrix for Compton scattering using 

The potential [V/ia^r] is calculated using the results of fl]], with only one exception: at the 
rjNN-vertex we now use pseudoscalar (PS) instead of pseudovector (PV) coupling. From 
eta photoproduction p9[ it is known that PS coupling leads to a much better description of 



the data. In purely hadronic reactions this PS PV difference is hardly visible because of 
the suppresion of the Born terms due to the hadronic form factor. But since, at least for the 
Ohta-prescription, the influence of the form factor is removed by restoring gauge invariance, 
in photoproduction the contribution of the Born terms is enhanced as compared to others. 
For kaon photoproduction the situation is not so clear [^], so we will use both PS and PV 
coupling at the KNA, S-vertices in the fits. 



D. Summary of the model 

Our model thus consists of a i^-matrix treatment of the following asymptotic channels: 
ttN, ttttN, riN, KA, and 7A^. The i^'-matrix elements are calculated from a Lagrangian given 
in for the hadronic couplings and in the present paper for the electromagnetic couplings. 
Possible intermediate states are the nucleon and all three-star nucleon resonances with spins 
1/2 and 3/2 up to an invariant mass of \/L9 GeV; thus no nucleon resonances appear in 
the final states of a Feynman diagram. As a consequence of the i^-matrix approximation 
all intermediate particles propagate only on-shell. 

With these ingredients all Feynman-amplitudes are calculated taking into account the 
S-, u- and t-channel diagrams. The latter are used to consistently generate the background 
amplitudes thus eliminating the need to introduce separate ad-hoc background parametriza- 
tions. The iteration of the amplitude V contained in Eq. then leads to a consistent 
description both of resonance decay widths and rescattering in a hadronic reaction. For 
example, in a reaction with asymptotic channels — > rjN a whole chain of rescatterings, 
like e.g., — > A^* — > ttN — > N** — > rjN , where N* and A^** are different nucleon reso- 
nances, is automatically generated. While the resonances in this way all acquire width, and 
are thus 'dressed', higher order contributions to the vertices are not taken into account and 
the extracted couplings therefore are 'bare' ones. 
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IV. RESULTS OF THE FITS 



The parameters of the model can be grouped into resonant and non-resonant ones. The 
first group is given by the masses and the hadronic and electromagnetic couplings of the 
resonances. Furthermore, it also contains the z-parameters of the spin-| resonances and the 
cutoffs A of the form factors at all vertices. In the case of the non-resonant background we 
have the couplings of the scalar and vector mesons to the nucleon and the cutoffs at those 
vertices. The vector meson decays are calculated with the couplings fixed to (|15). Only 
(7^^pO and g^riuj have been allowed to vary in the fits because the corresponding decay widths 
are only known to about 25%. 

As has been explained in Sec. |IIB| , only resonances with spin | and | have been taken 
into account. These are: 5ii(1535), 5ii(1650), Pii(1440), Pii(1710), ^13(1520), ^13(1700), 
^31(1620), P33(1232), P33(1600) and ^33(1700). Furthermore, the potential contains t- 
channel contributions from the hadronic and electromagnetic decays of the vector mesons 
p, uj, Qq and K*. For the coupling of the u to the nucleon, the values 

9U.NN = 7.98, K^NN = -0.12 (25) 
have been used P|J22]. The dependence of the electromagnetic couplings of the P33(1232) on 



the cuA^A^-couplings have been investigated in [|T0|,|T3|,|14[ . Mainly due to the cu-exchange con- 
tributions to the M]^{.^-multipole both couplings gi and g2 of the P33(1232) are infiuenced by 
variations in g^jMN and KujNN- Therefore, also the extracted E /M-iatio is somewhat sensitive 
to the values of those couplings. In principle one could go ahead and also determine g^jNN 
and Ki^NN through fits to the data. However, we have not chosen to do so, because we feel 
that the pion photoproduction data alone do not offer enough sensitivity to reliably deter- 
mine the cuA^A^-couplings. In view of the aim of this study to find a simultaneous description 
of all included channels the so induced error in the P33(1232)-couplings is acceptable. 



A. Reaction channels and database 

The hadronic database has been described in and consists of PWA results for both 
ttN ttN, ttttN [^^0 and cross section and polarization data for n^p rjn, K^K. Because 
of the much larger database used in the SM95-PWA as compared to the KA84-solution, we 
only use the SM95 data and the corresponding parameter set of |l[] in this calculation. 

• IP — ^ IP- Differential cross section data from various measurements 0] have been 



used in the fits. Furthermore, we include the LEGS data on photon polarizations 
the old data on the recoil polarization from Wada et al. |^T| have not been used here. 



Since, from the helicity couplings given in |T^, we expect sizable contributions from 
spin-| resonances (Di5(1675) and Pi5(1680)) not taken into account here, we fit the 
Compton data only for energies < 1.6 GeV. Only in this energy range we can be sure 
that all resonance contributions are taken into account. 

— > TzN: In this channel we use the single-energy data of the multipole analysis 
SP97 0. In principle, also another analysis (MA97 [Q) is available, based mainly 
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on the latest MAMI and ELSA data. This analysis, however, only covers the energy 
region < 1.35 GeV. Because of this restriction, we do not use it in our fits. Nevertheless, 
differences between the two solutions SP97 and MA97 and our fits are shown in some 
of the plots. In a later stage of the investigation it would be very interesting to perform 
restricted fits using the MA97 data. This could allow to investigate the dependence of 
the P33(1232) couplings on the multipole analysis used. 

Unfortunately, the spreading of the single-energy data SP97 is much larger than for 
the ttA^- analysis SM95. In order to further restrict the parameters, we also use the 
speeds calculated from the energy-dependent solution SP97 in our fits. However, some 
precautions are necessary when incorporating these data, in particular near the reso- 
nance positions of the Pii(1710) and the Di3(1700). Even though a resonant structure 
cannot be ruled out from the SP97 single-energy solution, the SP97 speed data are 
smooth in the vicinity of both resonances. In our calculation, however, we always find 
resonant structures, even for resonances coupling only weakly to 'jN. For energies 
near these resonances, therefore, no fit to the SP97 speeds is possible. For this reason 
the speed data for the resonant multipoles were not used near the Pii(1710) and the 
P'i3(1700) in our fit. 

7 A/" — > rjN: In the energy range below 1.54 GeV we only use the very precise data of 
Krusche et al. for the differential and total cross section. For energies above that 
only sparse data from different groups are available For the total cross section 
there are also data from ELSA |^ for eta electroproduction at very small fc^ (= 
—0.056 GeV^), but no differential cross sections. From measurements on deuterium 
targets also neutron to proton ratios have been extracted p^. Furthermore, a few 



target asymmetry data are available [^. The latest measurements at GRAAL on 



photon asymmetries, however, are not yet published. 

7P — > K~^A: Here, the best data come from the SAPHIR experiment p9. The older 



measurements of differential cross sections and A-polarizations have been carefully 



investigated by Adelseck and Saghai [^]. Because of systematic deviations of certain 
datasets, the errorbars of these data have been enlarged. In our fits we also use these 
newly assigned errors. 

The fits are labeled (as in |l[]) by the vriV-PWA used to determine the hadronic parameters 
and the type of form factors for the s- and t-channel resonances. An additional number 
indicates the method used to gauge the Born contributions: 1 - Ohta's method with 



fixed hadronic parameters, 2 - Ohta's method with all parameters fitted and 3 - Haberzettl's 



method with fixed a^'s (cf. Eq. (p3|)). Furthermore, from the fits performed for the 
hadronic channels, it is obvious that the exponential form Fg leads to larger values of as 
compared to the other form factors Fp and Ff. Therefore, we do not use the parameterization 
Fe for the case of photon-induced reactions. Finally, since the fits performed in [Q using Fp 
and Ft lead to very similar descriptions of the data, we limit ourselves to fits starting from 
the parameter set SM95-pt (given in |l| and also in the Tables |ITT| - 0). Here one must keep 



in mind that the use of different parameterizations for the form factors introduces a source 
of systematic error that can be of comparable size to the statistical uncertainty induced by 
the error of the data. 
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In summary, the fit results are then labeled by SM95-pt-l, SM95-pt-2 and SM95-pt-3, 
standing for the SM95 partial wave analysis of vrA^ scattering of the Virginia group 0, all 
using the form factor Fp from Eq. for all vertices with propagating hadrons and Ft from 
the same equation for the t-channel diagrams. They employ Ohta's gauge fixing method 
with hadronic parameters determined from a fit to hadronic channels alone, Ohta's gauge 
fixing method with all parameters refitted to all channels, including the photonuclear ones, 
and Haberzettl's gauge fixing method, respectively. 



B. Fit with fixed hadronic parameters 

In a first fit, we allowed only the electromagnetic couplings to vary. All other values 
of masses and decay widths were taken from the parameter set SM95-pt of Using the 
procedure of Ohta to restore gauge invariance, the parameters have been determined by a 
simultaneous fit to the full database as described above. 

In Figs, m - ^ we show the results of this fit as dotted lines, together with the other fits 
described below. From the plots it is clear that the experimental data in all channels can 



be reproduced rather well. The improvement over our old non-unitary calculation [|22| using 
the T-matrix approximation (also shown in Figs. [T| - ^ is obvious. 

a. Comparison to a T-matrix calculation One of the most noticeable differences between 



the calculations performed here and those in Ref. ||2^ is the improvement in the description 
of the Ei^^- and M^_-multipoles. In the case of it is well known that only a correct 
treatment of the rescattering allows a quantitative description of this channel. The reason 
for this can best be seen in Fig. where a calculation with both P33(1232)A^7-couplings set 
to zero is shown. Even with no direct coupling to the resonance, the structure of the data in 
the £'i{.^-multipole can already be reproduced quite well. This shows that the rescattering 
is responsible for the shape of this multipole and not the direct excitation of the A. From 



this it is obvious, why our old model p2[ failed in describing this multipole 



In we have speculated that the same might be true for the M|'_ -multipole, but 
Fig. |10| shows that this is not the case. The direct coupling of the resonance is essential 
to describe the data for both -Eg-" ^^"^ M^_-multipoles. Therefore, the poor fit in the 
old calculation was obviously driven by the contributions of the Z)i3(1520) to some other 
multipole. These were most likely the offshell contributions that are not treated correctly 
in the T-matrix approximation. One important deficiency in this approximation is the 
appearance of spurious resonance-like structures (e.g. -Ef^, Mf^, E^,^ and Mfi^). These 
are induced by the offshell contributions of the spin-| resonances. As has been demonstrated 
in , these structures are an artifact of the T-matrix approximation and do not appear in a 
unitary calculation. Therefore, the investigation of offshell contributions of spin-| resonances 
and the corresponding z-parameters, as has been done recently by Mizutani et al. is 
not very meaningful in a T-matrix calculation. Without dynamical rescattering, the z- 
parameters are mainly adjusted to minimize the induced off shell-structures and reveal only 
little about the nature of the spin-| resonances. From what we have said about the M^_- 
multipole, it furthermore seems that also some resonance couplings are influenced by this 
effect. 

Also, the dynamical generation of the imaginary parts of the amplitudes leads to an 
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improved fit. Especially in cases where the Born terms dominate the amplitude, the old 
calculation did not generate the correct imaginary part {Eq^, £'f^ and Mf^). This is also 
easily understood, since in the T-matrix approach the Born terms are purely real. 

b. Compton scattering This is the one of the first attempts to calculate Compton 
scattering in a dynamical model beyond the A resonance. Therefore, we are for the first 
time able to check if the data on Compton scattering and meson photoproduction can be 
described using the same helicity couplings for the various resonances. As can be seen 
from Fig. we are able to fit the available data on Compton scattering very well. Both 
the differential cross section and the photon asymmetry S are reproduced over the whole 
energy range. From this we conclude that the Compton data are indeed compatible to the 
experimental results for the photoproduction channels. 

The main contributions to the cross section come from the Born s- and w-channel dia- 



grams and the resonances P33(1232) and Di3(1520). This can be seen from Fig. |TT], where a 
decomposition of the differential cross section into the individual contributions is shown. It 
is also obvious that the vr^ and t] t-channel diagrams have a small influence under backward 
angles only. For energies below 1.6 GeV all other resonance contributions could be safely 
neglected, none of them exceeds 5 nb/sr. 

Furthermore, in our calculation there is no need for an additional attenuation factor for 
the Born terms, as introduced by Ishii et al. pl| {x = cos 9): 



^Born — ^Born^ '^^^ ^\ (26) 

with a free parameter C fitted to the Compton data. The strong backward peaking of the 
Born contributions is an artifact of the T-matrix approximation employed by Ishii et al. and 
does not persist once the amplitudes are properly unitarized. To illustrate the difference 



between the K- and T-matrix calculation, we show in Fig. |T2| the Born s- and -u-channel 
contribution to the differential cross section employing both K- and T-matrix approximation. 
The inclusion of mainly 7iN rescattering leads to an enhancement of the cross section under 
forward angles and to the abovementioned reduction under backward angles. Therefore, we 
are able to fit the Compton data without an additional factor e~'"^^~^\ 

For the photon asymmetry S we also show the results of the isobar model of Wada et al. 



31| in the lower part of Fig. ^. Obviously, this ansatz is not able to describe the data. Both 
magnitude and shape are in disagreement with the experimental results. The dispersion 
relation calculations of LVov on the other hand, can reproduce the polarization data 
very well. In this region they practically coincide with our results, whereas for energies y/s ~ 
1.08 GeV both approaches differ by a factor of two. 

This observation allows us to investigate the validity of the i^-matrix approximation 
in some detail. The main difference to the dispersion relation calculation performed in 



1 4^ is the onshell-approximation for Gbs in Eq- (0)- Therefore, in our calculation the 
7tN intermediate state does not contribute below rriN + m-j^ ~ 1.08 GeV. Taking also the 
offshell-propagation into account (as the dispersion relation implicitly does), the Compton 
scattering 'sees' the vrA^-channel already for lower energies. This is the main reason for the 
different results in S within about 50 MeV around the vrA^-threshold. We thus conclude 
that the offshell 'tails' of the propagator Gbs do not extend much further than ^ 50 MeV. 
This observation is in agreement with the results of Pearce and Jennings [H], who found that 
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a rather soft offshell-cutoff in Gbs is needed (A ^ 300 MeV) to describe the TrA^-phasehifts. 
Thus, it seems that the ii"-matrix approximation yields rehable resuhs for energies not too 
close to a meson-production threshold. 

c. Eta photoproduction Looking at the overall result from our fit, we find that all major 
structures of the data are also visible in the calculation. Only for a few channels significant 
deviations from the data can be seen. The most prominent of these can be found in eta 
photoproduction below 1.6 GeV (cf. Figs. ^and|^). As we have pointed out earlier, especially 
in this region high precision measurements are already available |3^. Since the forthcoming 
experiments should yield data with comparable quality, the eta photoproduction can be seen 
as a 'testing ground' for all models that try to describe photon-induced reactions. Only the 
full description of this data in all details might allow the unambiguous extraction of the 
5*11(1535) resonance parameters. 

From the differential cross section it is clear that mainly the absolute magnitude is too 
small for energies below 1.5 GeV, whereas the isotropy is well reproduced. Therefore, an 
increase of the S'ii(1535)p7-coupling alone would not improve the overall fit, since it would 
lead to a drastic over-prediction of the data for energies above 1.5 GeV. From this we 
conclude that a change in the energy dependence of the resonance contribution is needed for 
a better fit in this channel. Such a change can only result from a variation of the hadronic 
masses and couplings; the coupling to the photon mainly influence the magnitude and not 
the shape of the 5*11(1535) contribution. 

This observation coincides with the fact that the poor 7c~p — > rjn data were responsible 
for the spread of the 5*11(1535) parameters between the different fits carried out in |jl|. Also 
the smaller ?7A^-couplings of the other resonances could not be extracted reliably. In the 
moment, it seems that the eta photoproduction imposes much stricter constraints on the 
resonance parameters, as the purely hadronic data does. This clearly shows that new precise 
photoproduction measurements need to be accompanied by improved hadronic data as well. 
Otherwise, extractions of resonance parameters will always be handicapped by the quality 
of the hadronic database. 

d. Kaon photoproduction In the other reaction channels this problem does not show up, 
mainly because of the lack of high precision data. Only in the case of kaon photoproduction 
for energies around 1.75 GeV we have indications of systematic deviations in backward 
directions. Here the cross section is dominated by the Born contributions, since Qkan is 
rather large (~ —6, e.g. compared to Qtinn ~ 1-2). In the hadronic channels only the 
product of coupling constant and hadronic form factor enters, which is much smaller (ps; 
—2.5 in the case of Qkan)- Additionally, n^p K^A is for higher energies dominated by 
the i^*-exchange in the t-channel. So the Born contribution, and thus Qkan, is not well 
determined by the hadronic data. 

In the kaon photoproduction, however, Qkan plays a dominant role, since in the Ohta 
prescription the influence of the hadronic form factor is completely removed. Furthermore, 
the contribution from the charge of the proton is not cancelled by a similar u-channel contri- 
bution, as it is the case in 77° production. Because of the hidden strangeness in the KA final 
state, we have a A (or E) propagating in the crossed diagram that only exhibits magnetic 
coupling (as in tt~^ production). Since all major contributions are therefore fixed from the 
outset, the fit could only be improved by reducing the photon coupling of the 5*ii(1650). The 
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resulting value (A^^^ = 31x10 ^ GeV ^/^) is significantly smaller than the number deduced 
from pion photoproduction = 69 x lO^^ GeV-^/^, 0). In our fit to the combined 

data of both channels the Eq_^_ data on pion photoproduction obviously do not play such 
an important role as the kaon photoproduction data because of the large uncertainty of the 
former in the region of the S'ii(1650). 

Already with fixed hadronic parameters we obtain a good overall fit. From the observed 
deviations it is clear that a further improvement can only be achieved by simultaneously 
varying some of the hadronic parameters as well. Before we show the results for such com- 
bined fits, we want to stress again that already with fixed hadronic parameters a reasonable 
description of all data is possible. Especially due to the dynamical rescattering, the main 
shortcomings of a T-matrix calculation have been resolved. 

C. Global fit using Oiita's prescription 

In this section we now discuss the results of a global fit to all hadronic and photon-induced 
channels in which also the hadronic parameters are allowed to vary. Ohta's prescription was 
used to gauge the hadronic form factors. The results are also shown in Figs. |] - 

a. Compton scattering Looking at the plots for the different reaction channels, we in 
general find only a slight improvement using SM95-pt-2. In the case of Compton scattering, 
the fit with fixed hadronic parameters already describes the data rather well, so that the 
new fit leads only to a relatively small decrease in x^-y (7.15 — > 5.20). Since the main 
contributions here come from the Born terms and the P33(1232) and Di3(1520) resonance, 
the changes in the differential cross sections can easily be explained by the slightly different 
helicity couplings found in both fits. For the P33(1232) both A1/2 and A3/2 are reduced and 
lead to the observed reduction for energies up to 1.3 GeV. In contrast to this, the increase in 
the Di3(1520) helicities increases the interference with the other contributions to Compton 
scattering. Therefore, the cross section is reduced slightly in this well. 

b. Pion photoproduction For the pion photoproduction the reduction of xl^^/ DF is due 

to a the better fit of the Mj'^.^-multipole. The increase in x^^^/DF comes mainly from the 
S'li-channel, since the S'ii(1535)-parameters exhibit the largest changes as compared to the 
SM95-pt-l values. Except from this we note only minor changes, mainly for channels where 
the background is dominated by the Born contribution (e.g. Mi_, and Accordingly, 
the values for the helicity couplings we extract are very similar for both fits SM95-pt-l and 
SM95-pt-2. In general, the agreement to the PDG-values is quite good. Serious discrepancies 
we find for the S'ii(1650) (for the reasons discussed in the last section) and for both the 
Pi3(1720) and D33(1700). That we find no agreement in the case of the Di3(1700) comes 
as no surprise, keeping in mind that this state is not well established and is found at rather 
different energies in the different analyses. Furthermore, the helicity couplings are known 
to be small and have very large errorbars. 

For both the Pi3(1720) and 1^33(1 700) the background is mainly due to the Born terms. 
As can be seen from Im(Pf_|_), Im(M"^) and Re(p2-^), this background is too large for higher 
energies. Accordingly, the helicity couplings of both resonances are adjusted to compensate 
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this contribution. Especially for the 1^33(1700) it is obvious that no good fit to the multipole 
data is possible with this large background. 

c. Eta photoproduction The drastic reduction of x^r^/DF (6.09 3.00) is accompanied 
by only a small increase of x^^/DF (1.73 1.95), so that we have an overall decrease for 
both channels (4.25 2.56). Furthermore, the dramatic increase of the Di3(1520)A^r7-decay 
width again shows the importance of a global fit to the full dataset. Fits to the hadronic 
data alone always yield very small values for this decay (~ 10 keV), whereas the combined 
fits are much closer [k, 50 keV) to the values found elsewhere (e.g. ~ 130 keV in |]18|), if 
one takes into account the lower mass of the Di3(1520) found here. In addition also the 
rjN -dec&y width of the 5*11(1535) significantly increases. As can be seen from Figs. ^ and |^, 
this increase is driven by the better fit to the cross sections close to threshold. 

Since we now have a reasonable agreement with the precise data of Krusche et al. |p4|| , 
we can turn to the extracted polarization observables. The results for the polarized photon 
asymmetry E, the recoil nucleon polarization V and the polarized target asymmetry T are 
shown in Fig. |], together with the few datapoints available and the calculations of Knochlein 
et al. for the photon asymmetries. The agreement with their calculation up to 1.6 GeV 
is obvious. Since S is dominated by the Di3(1520) contribution it seems that the r]N- 
coupling of this resonance is already well determined by the differential cross section. In 
contrast to this, we are not able to reproduce the target asymmetries for the low energies. 
At threshold none of our fits shows the measured forward-backward asymmetry. This is in 
agreement with the results of [^] that practically coincide with ours in this energy region. 
In a recent analysis Tiator and Knochlein have investigated the target asymmetry in 
a more phenomenological approach and have shown that a reasonable description of all 
data is only possible if one assumes a rather large, energy-dependent phase between the 
5*11- and Di3-contributions; such a phase is obviously not present in our results. For the 
higher energies we find no consistent results for all three polarization observables. Above 
the 5'ii(1535)-resonance the small contributions of the various resonances coupling to rjN 
interfere strongly with each other and with the background contributions. Further detailed 
investigations have to show, if the sensitivity of S, P and T to small contributions can be 
used to uniquely disentangle the different resonances. In the moment we can only observe 
that the polarization observables are not well determined from the fit to the differential cross 
section alone. 

In addition to S, P and T we also show our results for the neutron/proton ratios of the 
differential cross sections. From Fig. |^ it can be seen that the few data points do not put 
strong constraints on the fits because of their large error bars. Nevertheless, it seems as 
if the helicity coupling of the neutral 5*11(1535) deduced from the i?Q^-multipole from pion 
photoproduction is too small to yield the measured neutron/proton ratios. 

We find large variations of dan/ dap for the higher energies. From the differential cross 
section it can be seen that dap/ dVt is rather small for forward and backward angles. There- 
fore, we are extremely sensitive to the exact numbers obtained in the calculation in this 
regions. This indicates that dan/ dap is not a good quantity to investigate if one of the cross 
sections is close to zero. For example, even if there would be data for both channels under 
9 = 180° with an accuracy comparable to the results of Krusche et al. (~ 0.01 fib/sr), we 
could still vary dan/ dap by an order of magnitude without loosing the fit to the differential 
cross sections. Therefore, we define an isospin asymmetry X similar to the polarizations: 
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J 




(27) 



da,n + dan ' 



which is hmited to — 1 < X < 1. Besides this more technical advantage it also has a simple 
interpretation in terms of isoscalar/isovector couplings, provided one contribution to the 
amplitudes is dominant: 



Obviously, I vanishes only if either Qs or vanish. If the coupling to the proton {qs + g^) 
or to the neutron [qs — g^) vanish, I takes on its maximum values =f1. Furthermore, in the 
case of one dominant amplitude X should be rather isotropic. 

The results for X are shown in Fig. ^ It can be clearly seen that two different pro- 
duction mechanisms for forward and backward angles develop above 1.6 GeV. Below this 
energy the amplitude is dominated by the S'ii(1535)-contribution. The positive value for 
the higher energies in forward directions can be understood from the then dominant p- and 
cu-contributions to eta photoproduction. Since both add up for the proton case and have 
opposite sign for the production on the neutron, we would clearly expect X > 0. Even the 
magnitude can be explained in this simple picture: taking g^rjpo = 1 and g-yj^^ = 0.3 we 
readily obtain X = 0.55 from Eq. Under backward angles the situation is more compli- 
cated, since we there do not have one dominant contribution to the cross section. Here the 
Born terms, which are determined by cn and km, play an important role. Already these two 
couplings alone have different decompositions into isoscalar and isovector couplings. The 
same is true for the small contributions of the nucleon resonances. Therefore, one would not 
expect a simple explanation for the extracted values of X in this region. 

In summary, we find that the quality of the fits naturally improves, once we allow the 
hadronic parameters to readjust. The improvement is most significant in the eta photopro- 
duction. Mainly the resonances S'ii(1535) and Di3(1520) are affected by such a readjustment. 
For the other hadronic parameters, the masses of the resonances and the branching fractions 
do not change in a global fit using Ohta's method (except for 5'ii(1535) and Di3(1520), as 
explained above). The partial decay widths vary, but only in some cases (P33(1600), F^jv 
346 MeV 494 MeV and ^33(1700), F^jv 477 MeV ^ 337 MeV) the changes are signif- 
icant. Furthermore, we have found that the fits cannot be improved in channels that are 
dominated by Born contributions. Since we have used Ohta's prescription to restore gauge 
invariance for both fits, the contributions of the Born terms are nearly unchanged. From 
this we conclude that a further improvement is only possible, if the couplings to the charge 
of the nucleon, the pion and kaon are also changed by the inclusion of a form factor. 



From the x^-values (Table |T|) it is obvious that already a fit with fixed a^'s leads to a 
further improvement. We find that mainly the better fit to the Compton scattering and the 
pion production data are responsible for this, whereas the x^-values for the other channels 
remain fairly constant. 



X 



{.gs + gvf + {gs -gvf 



(28) 



D. Global fit using Haberzettl's prescription 
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a. Compton scattering The main improvement can be found in the fits to the differential 
cross sections at higher energies (cf. Fig. |^). Looking at the photon asymmetry S, one is 
tempted to conclude that the fit SM95-pt-3 is worse than SM95-pt-l,2. The solution to this 
puzzle is that, in the x^-analysis, the S data around 1.2 GeV are more important than the 
other points, since their errorbars are much smaller. Since these three points are reproduced 
better in SM95-pt-3, the x^-y does not increase, even though the slope of the data seems to 
be described better using the other two fits. 

b. Pion photoproduction In the pion photoproduction the most significant changes can 
be seen in the i?i+,Mi+- and i?2-,M2_-mutipoles. This can be understood from the ansatz 

= l/3(-F(s) + F{u) + F{t)) (pBI). The Born contributions to the s-wave multipoles i?o+ 
for example are mainly affected by F{s), since F{u) and F{t) induce angular-dependent 
modifications. Therefore, the changes in the s-wave contributions are not very large for 
energies < 1.5 GeV. Also the larger changes in the helicity amplitudes of the P/3- and Dj^ 
resonances (cf. Table [VID show that the Sn- and P/i-channels are affected only for higher 
energies (e.g. for the Pii(1710)). 

An interesting effect can be see for the changes in the -E'2-^- and M2/^^-multipoles as 
compared to M^_. For the first two we have already concluded in Sec. [IV C] that an improved 
description might only be found by changing the Born contributions. This can now be 
confirmed using Haberzettl's method. Also, the helicity couplings of the 1^33(1700) are now 
in somewhat better agreement with the PDG- values, as can be seen in Table [VI I[ Obviously, 

led to a significant reduction of the non-resonant background (cf. Fig. |]). Since P^- does 
not depend on isospin, we expect to have a similar reduction for £"2!" and Mfi". That this 
is indeed the case can be seen in Figs. || and H. In all four multipoles the agreement to the 
data for energies < 1.5 GeV is reduced due to the smaller background. The readjustment 
of the P>i3(1520)-parameters in fit SM95-pt-3 then results in some deviations, mainly in the 
M^_-multipole, since there the errorbars are largest. From this we conclude that SM95-pt-3 
represents a compromise between the improvement in the D33(1700)-case and the larger 
deviations for the multipoles containing the Di3(1520). 

Additionally, in all three fits we find that we overestimate the -Eg^i^-multipoles for energies 
around 1.3 GeV (Figs. |^ and^). Only part of this is due to the S'ii(1535), as can be seen from 
Fig. |T3|, where the results for Re(i?o+) shown with and without this resonance. From this 
it seems that the background is too large in this energy region. It is interesting to note that 
a similar discrepancy was also found in the i^'-matrix calculation of Deutsch-Sauermann et 
al. [jl2|, where an even larger non- resonant contribution was found (cf. Fig. |l^). Since the 



background for this energies is dominated by the Born terms, it might be possible to find a 
better description of the data if the parameters in (p3[ ) were also allowed to vary. 

c. Eta photoproduction Since QrjNN is small compared to Qt^nn, we find only minor 
changes in the case of eta photoproduction. The differential cross section in backward 
directions is larger for energies < 1.7 GeV; above this energy we have a reduction as compared 
to SM95-pt-l,2. This is due to the changes in the interference of the Born terms with the 
p- and C(j-contributions. Since we have no dominant resonance in this energy range, these 
changes can be easily observed. 

d. Kaon photoproduction For the kaon channels, both x^k (3.91 — > 4.09) and xix 
(3.77 —>■ 4.21) increase slightly as compared to SM95-pt-2. Obviously, in the case of kaon 
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photoproduction the effect of Fx is largely compensated by the increase of Qkna (—6.25 
—8.65, cf. Table |IT|). Since Fx = l/2{F{s) +F{t)), this compensation cannot be complete. 
Because of the t-dependence of F^ a small reduction under backward angles and a similar 
increase in forward directions is expected. The calculated differential cross sections indeed 
show this behavior (Fig. P). It can be seen that the use of Fx does not lead to an improved 
description of the cross sections. 

In contrast to this, the A-polarizations can clearly be reproduced better using SM95-pt-3. 
Especially close to threshold we find the polarization to have the right sign and magnitude, 
in contrast to the other fits SM95-pt-l,2. Responsible for this improvement are not the 
changes in either Born or S'ii(1650)-parameters, but mainly the Pii(1710)p7-coupling that 
increased by a factor of two in SM95-pt-3 (cf. Table [VI| ). 



In summary, we find that Haberzettl's method, aside from its theoretical appeal in that it 
incorporates the physical constraints on all momenta, also leads to a better fit, mainly in the 
Compton scattering and photo-pion channels. In a calculation, in which also the parameters 



tti in eq. are allowed to vary, even more significant improvements can be obtained p? 



V. PARAMETERS AND COUPLINGS 

After the more phenomenological discussion in the last section we now want to focus on 
the extracted parameters. To this end we first investigate the non-resonant couplings and 
after that the resonance parameters as found in the different fits. 



A. Background parameters 



In Table |ITl| we list the final values for the couplings of the mesons to the N and A and 
also those for g^^ipo and g-yrjoj. The other photon-decay couplings of the mesons have been 
kept at the values given in ([T3|). For these the errors deduced from the uncertainties in the 
decay width are of the order ~ 5%, whereas for g^npO,oj we have ~ 25%. Furthermore, only 
in the eta photoproduction we have some sensitivity on the background couplings also at 
higher energies because there is no dominant resonance contribution. From previous studies 
29| , pn]jr^ it is also know that in this channel we have a large cancellation between the Born 



contributions and the p and uj t-channel exchanges. This also enhances the sensitivity of the 
fits onto the parameters g-yripO^^. 

For gj^NN and gr,NN obviously all fits yield very similar results, comparable to the SM95- 
pt-1 values, which have been deduced in |^ from the hadronic data alone. For g-i^NN this 
comes as no surprise, since this value has already been extracted many times consistently 
from hadronic and photon-induced reactions. 

In the case of grjNN other groups find somewhat larger values than the ones deduced 
here: 2.24 from eta photoproduction p9|, 6-9 from NN potentials In our analysis 

the main sensitivity comes from the data under backward angles for both n^p rjn and 
7P — > rip. Interestingly, the values we find are even smaller than those from other fits to 
eta photoproduction: Benmerrouche et al. find gr,NN ~ 5 in a T-matrix calculation using 
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effective Lagrangians |TT|. In contrast to that Tiator et al. p9| deduce 2.24 from a model 



that tries to incorporate unitarity by the use of energy dependent phases at the RN^f- 
vertices. Since we also find somewhat smaller cross sections (~ 10 %) using the T-matrix 
approximation, our even smaller values grjNN ~ 1-0 might be due to rescattering effects. 
Especially, since we have the above mentioned cancellations between different non-resonant 
contributions, the extracted values for the background couplings are rather sensitive to the 
approximation used. 

For the Qkna we have a totally different situation. Here, the fits using Haberzettl's 
prescription find a much larger value than those using no form factor at the charge contri- 
butions. However, as we have already pointed out in Sec. [IV D| , the effective coupling g ■ F 
is of the same order for all fits, e.g. g ■ F = —3.79 for SM95-pt-3. This coupling is therefore 
a very good example that it does not make sense to compare the bare couplings deduced 
from different ansatzes/reactions/models once form factors have been taken into account. 

Since we are now able to describe both the meson- and photon-induced production of 
KA using rather small values gxNA (~ —6), it seems that the discrepancy to the SU(3) 
predictions (~ —(10.3 - 16.7)) cannot be removed easily. But as can be seen from SM95- 
pt-3, only the effective coupling is determined by the fits. Therefore, the SU(3)-values can 
of course be used, as long as one introduces a suitable form factor F. This is clearly not 
very satisfactory, since it renders the whole procedure of determining gxNA using SU(3) 
questionable. Only from a microscopic model for the form factor one could judge, if a bare 
coupling gxNA compatible with SU(3) would lead to reasonable fits to the data. 

One other solution to this problem was sought in the use of PS- instead of PV-coupling 
at the i^A^A-vertex. The investigation of kaon photoproduction in non-unitary models (see 
30| for a detailed discussion) suggested that PS-coupling leads to larger values for gxNA, in 



better agreement with SU(3). To check whether this conclusion also holds in a multichannel 
calculation, we performed a fit starting from SM95-pt-2 employing PS-coupling. We only 
show the most important results of this fit in Table 11^. From the x^-values we can see that 
both PS and PV yield fits of similar quality. However, we do not find a significant increase of 
gxNA (—6.8 instead of —6.2). Already the 7r~p —>■ K^A data alone are not compatible with 
\gKNA\ > 10, even though the contribution of the Born terms is suppressed by the hadronic 
form factors 



B. Resonance parameters 

Finally, we discuss the parameters of the nucleon resonances as found in the fits. These 
are collected in Tables |V| - |V11| , where also the PDG- values and the results of Arndt et al. 
for the various helicity couplings are listed. We do not show the results of other models 
for the purely hadronic parameters. These, and a detailed discussion, can be found in [|l|. 
As can be seen from Tables |VT| - |V14 the agreement of the helicity couplings deduced in this 
work to the PDG-values and the values given by Arndt et al. is quite reasonable. The most 
important deviations {Ay2 of the 511(1535), 5*11(1650) and the Di3(1520) and also for A1/2 
of the 1)33(1700)) have already been addressed in the previous sections. They have been 
related to the use of additional data besides the pion photoproduction multipoles and to the 
influence of the gauge prescription used in the fits. 
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/Sii: For the 5*11(1535) all three fits lead to helicity couplings much larger than the 
PDG-values. A similar discrepancy was also found in the extraction of A^^^ from eta pho- 
toproduction |^6|. There resonance parameters very similar to ours have been found: 

niR = 1.544 GeV 
Ttot = 212 ± 20 MeV 

= (120 ± 26) X 10-3 GeV-^/^ (29) 

Furthermore, it has been demonstrated that larger values for mji also lead to a larger width. 
This trend can clearly be found in our results as well. We therefore confirm the findings of 
| 46| that the eta photoproduction data can only be explained using helicity couplings larger 



than the PDG-value. 

Since the PDG-value has been extracted using pion photoproduction data alone, one 
might think of an inconsistency between both datasets. Our results show that this is not 
the case and that both reactions can be described using a large Ay 2- This observation was 
first made by Deutsch-Sauermann et al. ||T^. There the value of Ay^ = 102 xlO"^ GeV~^/^ 



was extracted by a combined fit to pion and eta photoproduction. For this conclusion the 
treatment of the rescattering seems to be important, since in a T-matrix calculation using 
effective Lagrangians a smaller helicity coupling was deduced (^4^/2 = 87x 10^^ GeV"^/^, 



2^ ) from the pion data. 

The values for A^^^ of the 5*11(1650) are found to be smaller than the PDG-values. It 
is mainly determined by the ii'A-channel, as has been discussed in the last section. Here 
indeed it seems that this small value is in contradiction to the pion photoproduction. Unfor- 
tunately, the data on the £'o+-multipole are not very good in this energy range. Therefore, 
they do not constraint A^^ very much. Especially a better determination of Im(£'o+) would 
help to clarify this situation. 



Pii: In the case of the Pii(1440) the fits SM95-pt-l,2 agree very well with the values 
obtained elsewhere. Only for SM95-pt-3 we find a somewhat larger coupling {Ay 2 = — 84x 
10^'^ GeV^^/^). Nevertheless, the fit to the Mf_-multipole is still as good as for the other 
two fits. The change in Ay2 is therefore be due to the use of Haberzettl's prescription, which 
leads to a reduction of the Born contributions. 

For the second resonance we also find larger couplings once we allow for a residual form 
factor F. Unfortunately, the data for pion photoproduction are not good enough to con- 
strain the fit. As we have discussed in Sec. |I V D| , the increase of Ay2 in SM95-pt-3 is mainly 
driven by the data on the A-polarization in kaon photoproduction. 

P13: In this channel we are not able to describe the imaginary part of the multipoles. 
There are indications for a second state with these quantum numbers (Pi3(1879), |^), which 
was not taken into account here. Maybe the inclusion of an additional resonance would lead 
to a better fit. Accordingly, the couplings that we find are not well determined and vary 
between the different fits. 

Dis'. For the Di3(1520) the most obvious deviation from the PDG-values is for Ay2 (-6 
in SM95-pt-l as compared to —24 ±9). In the multipoles this shows up in the Mf.-channel, 
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where we miss the imaginary part by roughly 20%. This disagreement is clearly driven by 
the simultaneous description of Compton scattering. This is illustrated in Fig. |14| where we 
show the Compton result using the PDG helicity couplings for the Di3(1520). With these 
values the Compton data are overestimated by a factor of 2-3, while with our couplings a 
good fit to the data is obtained. Also shown in Fig. |T^ is the cross section using only the 
P33(1232). It can be seen that this resonance gives a divergent contribution for energies > 
1.8 GeV. This is due to the rather large cutoff ~ 4.0 GeV), see Table |V111| ) obtained in 
all fits, which use the Compton data only up to 1.6 GeV. As a consequence, the P33(1232) 
gives a large background already at lower energies, accounting for ~ 50% of the cross section 
in the Di3(1520)-region. In this energy range the s- and w-channel diagram of the P33(1232) 
lead to a comparable cross section. Therefore, the small value of might be forced by 
the interference with an unreasonably large background. 

To investigate this question in more detail, we have performed a fit starting from SM95- 
pt-2 with a fixed cutoff Ay 2 = 1-1 GeV; the Compton scattering cross section obtained 

with this value is also shown in Fig. |1^. Mainly because of the data in the pion Mi[. - 
multipole the total increases to 8.70 and the values for the Compton channel and the 
pion photoproduction are found to be = 12.44, x^n = 10.88 (as compared to x^-y = 3.40, 
X^TT = 6.69 in SM95-pt-3). For the 1^13(1520) the resulting value of Ay2 changed from 3 to 
—6 and is therefore not much closer to the PDG-value —24 ± 9. 



Interestingly, in their old isobar model analysis of Compton scattering, Wada et al. ||31 
found A^^2 close to our values, namely A^^2 = —12.1x10^^ GeV^^/^. Thus it seems that 
this value is not dependent on the details of the model, but is forced by the Compton data; 
here clearly further investigations are necessary. 



The changes in the couplings for fit SM95-pt-3 have already been discussed in Sec. IV D 



There it was shown that the use of a residual form factor at the Born terms results in large 
changes of the background amplitudes. Especially for no good fit can be found in this 
case. 

As in the purely hadronic fits, we find the second Z}i3-resonance at rather high energies 
> 1.9 GeV. Since we limited our fits to this energy, we cannot meaningfully extract the 
parameters of this resonance. Therefore, the values given in Tables ^ and [V^ should only 
be seen as an indication that a second resonance exists in this energy range. 



S31: Even without a direct photon coupling of the 5'3i(1620) we can describe the mul- 
tipole data quite reasonably. As in the hadronic case, the background in this channel is 
dominated by the Born terms and the P33(1232) contribution. Therefore, the resonance pa- 
rameters of the 5*31(1620) are very sensitive to the z-parameters of the P33(1232). This can 
be seen from fit SM95-pt-3, where the sign change in ^41/2 is mainly driven by the change 
of Z2. 

P31: We do not include a resonance in this channel. Nevertheless, we are able to re- 
produce the Mj'^^-multipole up to energies ~ 1.7 GeV. Above that, in the imaginary part 
clearly the contribution of a higher lying state becomes visible. Therefore, we use the data 
only up to this energy. 
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Ps3- In the global fits the hadronic parameters P33(1232) tend to decrease even further. 
The numbers found are all at the lower end of the allowed region. An unconstrained fit would 
lead to even smaller values (m^j = 1.226 GeV, F^^r = 105 MeV), but the improvement in 

would be minimal. Since these small numbers can be understood in terms of the p 
contribution to ttN scattering at higher energies (as has been discussed in |1[]) we have 
chosen to limit the parameter range for the P33(1232). The numbers in Table [TVl, which are 
all at the lower bounds, therefore indicate that the hadronic parameters of the P33(1232) 
are still too small, even in a global fit to all reactions. 

All three fits yield somewhat smaller electromagnetic couplings for the P33(1232) than 
the PDG- values. As in the case of the Di3(1520), this is due to the inclusion of the Compton 
data. A fit without this channel would lead to somewhat larger couplings. Nevertheless, 
these changes are small. In this energy range both reactions can therefore be described by a 
single set of parameters that is in agreement with the values deduced from photoproduction 
alone. 

The £'/M-ratio deduced in our fits is 0: 

= -2.1(3)%, (30) 

where the number in brackets denotes the error in the last digit. There has been some debate 
about the exact value of this quantity, but it has been demonstrated that the different values 



extracted were due to the use of different datasets (see |^ for details). Using the SM97- 
PWA, Arndt et al. deduced Re/m = -1.5(5) %, the use of the MA97-PWA led to Re/m = 
—2.4(4) %. Wilbois et al. have pointed out that the E'/M-ratio extracted from dynamical 
models depends on the unitarization method used From the different ansatzes they 

found Re/m = —(0.7 - 5.7) %. In order to have a model-independent quantity it was 
proposed |^ to use a speed plot analysis to determine the i?/M-ratio at the pole of the 
P33(1232) (for a detailed discussion see PH]). From this one finds [ff^: 



SP97 : Re/m = -0.034(5) - 0.055(5)1, , . 

MA97 : Re/m = -0.035 - 0.046i. ^ ^ 

Using the same technique here for the numerical results of our fits, we extract an E / M-ratio 
of: 

i?^/^,^ = -0.022(3) - 0.002(2)1. (32) 

The obvious disagreement to the values given in (^1]) can be traced back to the imaginary 
part of the i?^{.^-multipole. As can be seen from Figs. ^ and |T0|, it does not rise steeply 
enough in our fits. 

For the second P33-resonance we find rather large values for ^3/2- Furthermore, in the 
fits SM95-pt-2,3 also the mass and width increase as compared to the purely hadronic fits. 
Both effects are driven by the fit to the M^^^-multipole. As can be seen from Fig. 0, for the 
higher energies only a change in the mass and width leads to a good description of the data. 
In the case only SM95-pt-3 yields a reasonable fit for energies > 1.5 GeV. However, 
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also in this channel the data are not good enough to extract the helicity couplings of the 
P33(1600) unambiguously. 



Dss'. As discussed in Sec. [IV D| , a satisfactory fit could only be found using a residual 
form factor Ft^. This again shows up in the extracted helicity couplings. Only for SM95-pt-3 
we find good agreement with the values obtained by Arndt et al. 0. 

2-parameters: We confirm the finding of that only the 2;-parameters of the P33(1232) 
and the Di3(1520) can be extracted reliably. Since these two resonances give large contri- 
butions to the background in the S- and P-waves, their magnitude can be determined 
independently form the resonances in these channels. 

Clearly visible in all cases is a dependence of the z^s on the gauge prescription used. 
This can be understood because the residual form factor F changes the non-resonant Born 
terms, so that the background determined by the 2;-parameters needs to readjust during the 
fits. Mainly Z2 is effected by this, whereas and Zi remain rather stable. The values for 
the P33(1232) are found to be: 

2^ = -(0.31-0.35), ;zi = -(0.27- 0.53), ^2 = -0.66 - 1.37, (33) 



and can be compared to the extractions by Davidson et al. |T0[ and by Olson and Osypowski 
| 50| who, however, do not use the second coupling of the P33(1232) to 7A^: 



Davidson et al.: Zt^ = —0.24, zi = —0.53, Z2 = 2.39, 

Olson and Osypowski: = -0.29(10), zi = 0.78(30). 

The values found here are obviously in good agreement with both extractions. 
For the Di3(1520) our results are given by: 



(34) 



z^= 0.31-0.35, z^= -(0.17-0.66), 2^= 0.57- 0.82, 

zi= -(0.26 - 0.32), Z2 = -(0.48 - 1.36). 



(35) 

Here no other systematic investigation is available. Fits to the eta photoproduction alone 



found no sensitivity to these parameters [11 . 



VI. SUMMARY AND CONCLUSIONS 

The aim of this paper was to extend the ii'-matrix calculation of to include photon- 
induced reactions. To this end we have introduced the final state into our model, leading 
to the new reaction channels jN — >■ jN, nN, rjN, KA. From a fit to the combined database 
of hadronic and photon-induced reactions three parametersets have been extracted. 

It was shown that keeping the hadronic parameters fixed to the values obtained in |jl| 
no satisfactory fit to the eta photoproduction data could be obtained. Therefore, at least 
in this channel, the photoproduction measurements are needed as additional input for the 
extraction of the resonance masses and width. 

The use of hadronic form factors at the 99 iVA^- vert ices violates gauge invariance so that 
counter terms have to be introduced to restore it. Unfortunately, these counter terms cannot 



26 



be constructed in an unambiguous way. To investigate this systematic uncertainty, we have 
performed fits using both Ohta's and Haberzettl's method. 

We find that only a few couphngs are sensitive to the gauge prescription used. Mainly 
the parameters determining the non-resonant contributions (e.g. the ^-parameters of the 
spin-| resonances) change between the fits. Furthermore, some partial waves are affected 
more than others. The reason for this is the additional angular dependence introduced by 
the residual form factor F{s,u,t) found in Haberzettl's ansatz. Even though we find it to 
be small, this model-dependence has to be kept in mind when extracting helicity couplings 
using effective Lagrangians. Only in a microscopic model for the hadronic form factors this 
problem might be resolved. 

We have shown that the Compton data provide an important additional source of infor- 
mation for the extraction of the electromagnetic coupling constant. The combined use of 
the Compton date and the pion photoproduction multipoles shows only one possible can- 
didate for an inconsistency between both data sets: the helicity couplings of the Di3(1520) 
are found to be smaller as compared to values extracted from pion multipoles alone. Here 
further investigations are necessary, since we also find a rather large background in the 
Compton amplitudes. For the P33(1232) a combined fit yields couplings ~ 5% smaller than 
the PDG-values. 

In the case of eta photoproduction the cross sections can be fitted rather well, but 
we are not able to reproduce the measured target asymmetry close to threshold. For the 
higher energies we find that especially the polarization observable are sensitive to the tails 
of the S'ii(1535) and Di3(1520). The interference of weekly coupling resonances with this 
background leads to large fluctuations between the different fits. This effect might be used 
to look for these resonances, since they are otherwise hard to detect from the differential 
cross section alone. 

In all fits we find a i^'A^A-coupling of half the size of the SU (3)- value. Since this results 
also persist using PS-coupling at the i^'A^ A- vertex, we conclude that with Qkna taken from 
SU(3) no fit to the combined hadronic and photoproduction data is possible. 

The basis for all these results is the fT-matrix method used in the present study. The 
main approximation in this method is that all intermediate state particles are put onshell; 
this violates causality since the connection between real and imaginary part of the propa- 
gator is lost. This approximation is clearly unphysical for bound states, but for scattering 
states it seems to be quite reliable. We conclude this from the cited studies of Pearce and 
Jennings [Q and Surya and Gross Further, more heuristically, the agreement of the 
resonance parameters extracted in the present study with those obtained from other uni- 
tary analyses (for smaller channel spaces and limited energy regimes) also shows that the 
ii'-matrix approximation is quite reliable. 

We have tried to obtain a more quantitative estimate for the validity of the i^'-matrix 
approximation from a comparison with dispersion relation calculations for Compton scat- 
tering. There we have found that both calculations agree very well except for energies ± 
50 MeV around the pion threshold. This is quite understandable in view of the discussion 
above: just above threshold the bound state behavior can still be felt. We thus expect the 
i^-matrix approximation to be resonable also for higher energies, except close to particle 
production thresholds. Any remaining shortcomings of this method also have to be seen 
in comparison to its inherent and practical advantages. First, because of the absence of 
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off-shell propagators no regularization of amplitudes is needed. There is therefore no need 
to renormalize, for example, coupling constants etc. Second, after a partial wave expansion 
of the scattering amplitude has been made, the original integral equation reduces to an 
algebraic equation which can more easily be solved. 

In the future it is clearly desirable to combine the multi-channel calculation performed 
here with approximations beyond the /^-matrix ansatz, as they have been developed e.g. by 
Surya and Gross and Sato and Lee []T1|. Work along these lines is currently under way. 



In summary, we were able to describe the combined dataset of hadronic and photon- 
induced reaction using the same set of parameters. This model can therefore be used for 
numerous detailed investigations that have not been possible before. 
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APPENDIX A: HADRONIC COUPLINGS 



In this section we list the hadronic couphngs used in our model. For a more detailed 



account, see [0]. Throughout this paper we use the notation of Bjorken and Drell [^. Four- 
momenta are denoted by a;, x is the absolute value of the corresponding three- momentum 
X. Furthermore, x is a. unit- vector in the direction of a;: x = x/x. 
For the nucleon, the following couplings have been used: 

^NR = - ^^^N'^^'^^{d^>^)N - gsNNs{NN) - gs^^s{ip*ip) 

- g.NNN {^i^v^ - ^v^V'"^ N - g,^^ [ip X {d,<f)] v^. (Al) 

Here ip denotes the asymptotic mesons vr, r] and a coupling to the C-meson is not taken 
into account, s and v are the intermediate scalar and vector mesons (oq, p and K*) and 
y^lu _ ^i^^M _ Qt^y^ is the field tensor of the vector mesons; is either a nucleon or a A 
spinor. For the /= 1-mesons (vr and p) if and need to be replaced by r • and t • m 
the 99, f A^A^-couplings and by (p and otherwise. 

For the Sn resonances we employ (pseudo-) scalar coupling to mesons and nucleons in 



accordance with |jT2[. The (pseudo-) vector coupling is used in the case of nucleon and all 



P/i resonances. This is done in order to reproduce the vrA^-scattering lengths without an 
additional sigma meson. For the spin-| resonances offshell parameters Zi are taken into 
account at all vertices and fitted to the available data. 
For the Sn and S'31 we therefore use: 

^NRy, = -Q^nrR T^N + h.c, (A2) 
and in the case of Pu and P31 the couplings are given by 

iil'R =t 11 I'M 

with the upper sign for positive parity. The vertex-operators F and F^ depend on the parity 
of the particles involved. For a meson with negative intrinsic parity coupling to two baryons 
with positive parity (e.g. ttNN) they are given by F = 175 and F^ = 757^, otherwise (e.g. 
7rA^S'ii(1535)) we have F = i and F^ = 7^. 

For the spin-| resonances the following coupling is used: 

Qaf.{z) = gaf. - ^(1 + 22;)7a7/., (A4) 

again with a vertex-operator F that is 1 for a particle with negative intrinsic parity and 75 
otherwise. 

For the isovector-mesons vr and (, (f in ( |A2| ) - ( [A4| ) needs to be replaced by r • <^ for 
/ = I resonances and by T • <^ otherwise. 
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The couplings constants can be derived from the decay widths using the following for- 
mulae (p denotes the three-momentum of the meson and nucleon, Ej^ and E^p the nucleon 
and meson energy, respectively): 

For spin- 1 resonances we have: 

PS-coupling : 

i ± — ibU — p r= 

PV-coupling : 

r ^Inr '^E^jE^E^ + p^) - mliE^ ± m^) 

L± — IbU — — — — p ^= . [AO) 

The upper sign corresponds to decays of resonances into mesons with opposite parity (e.g. 
Pii(1440) — >• TiN), the lower sign holds if both have the same parity (e.g. 5*11(1535) — > ttN). 
ISO is the isospin factor, it is equal to 3 for decays into mesons with isospin one, 1 otherwise. 
Spin- 1 resonances: 

Again, the upper sign is used if resonance and meson are of opposite parity. 



APPENDIX B: EXTRACTION OF PARTIAL WAVES AND MULTIPOLES 

In general, the extraction of the partial waves in hadronic reactions and the multipoles in 
the photon-induced channels is done similarly. Starting point is always the invariant matrix 
element M. fi and its connection to the scattering amplitude: 

Mfi = uip', s')du{p, s) = i^x^^Xi- (Bl) 

Here, O denotes the transition-operator for a given reaction, m and m' are the masses 
of the initial and final baryon, respectively. The decomposition of T in terms of partial 
waves/multipoles is usually straightforward, but tedious. Once the relations between the 
expansions of O and !F have been established, any contribution to the invariant matrix 
element can therefore be decomposed into partial waves/multipoles. In the following we 
want to list the relations O — > ^ for the three types of reactions that are described in our 
model. 

1. Meson nucleon scattering 

In this case, the invariant matrix elements A^/i and the scattering amplitudes JF are 
given by: 
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Mfi = uip, s')T {A + BQ) u{p, s) 
J^n = (Act -p' + B a -P), 



(B2) 



with r = 1 for mesons with equal parity in the intital and final state, F = 75 otherwise. 
Fe^n are the scattering amplitudes for the case of equal/non-equal parity, respectively. The 
relation between the amplitudes A, B and their counterparts A, B can easily be established 
P|JTB|, and in turn can be used to calculate the partial waves Ti± from the invariant matrix 
element (x = cos 6): 



A 
B 



{E' ±m'){E + m) 



(A + 5(v^-m)) 



{E' Tm'){E -m) 



Sttv^ 



(A-fi(v^ + m)), 



^qq' 



dx (APiix) + BPi^ix)). 



(B3) 
(B4) 



the upper sign holds for mesons with equal parity, m = (m' ± m)/2. The reader is 



In (I 

referred to 11161 for a detailed derivation of this relation. 



2. Photoproduction 



In this section we list the formulae needed for extracting the multipole amplitudes from 
the invariant matrix elements. This is by no means a complete listing of all the relations 



between the different set of amplitudes; a detailed account of this can be found in [43 



For the multipole decomposition of the meson photoproduction, we first need the relation 
between the invariant matrix elements A^/i and the CGLN invariants jFj ||52|. To this end 
we expand the Feynman amplitude as 



where 



iM fi = u{p\ s') AjMju{p, s) 
3=1 



Ml = -F^l^, 

M2 = 2T{{ep){kp')-{kp){ep')), 
Ms = T{^{kp)-Usp)), 
M, = T{^{kp')-Uep')). 



(B5) 



(B6) 



F = 75 for the production of mesons with negative parity, F = 1 otherwise. This decom- 
position differs from the one usually used in pion and eta photoproduction |Tl|^j52 



to 



simplify the inclusion of kaon photoproduction. The relation to the usual CGLN invariants 
Mi is given by: 
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Ml = Ml, M2 = -M2, M3 = -(M3 - M 



4;, 



M4 



-(M3 + M4) - 2m^Mi. (B7) 



In terms of Pauli spinors the scattering amplitude can be expressed using the jFj's [112 



jFp = icr • q cr • eT\ + cr-(fc X £)JF2 + icr • g cr • fc e • qTz + ie • 9^4, 
Tm = icT • + (T • q (T-{k X s)jF2 + icr • fc £ • qjF3 + icr • q e • qJ^^. 



(B8) 



Here, refer to the parity of the meson that is produced, k and q denote the three- 
momentum of the photon and the meson, respectively, i is a unit-vector in the direction of 

X. 

For our choice of amplitudes Mj, the relations between the CGLN invariants JFj and the 
A/s of Eq. 



are the following |^ 
^1 



^2 
^3 

^4 



k 

k 

47r' 



IE' ± m' 



2Vi 



IE' T m' 



-Ai-^A3- 



{kp') 



kq E'±m' 
4^ V 2Vi 



kq / =F m' |- 
4^V 2v/i 



2 ' + m 

s — m)A3 — A4 

{^/s + m)A3 + A^^ 



(B9) 



where the upper sign holds for mesons with negative parity, k and q are the absolute values 
of the three-momenta k and q. 

In order to unify the formulae for mesons of positive/negative parity, we introduce a 
notation different to the usual one, by labeling the multipoles using initial variables instead 
of final variables [^]. In Table | we list the relation between both notations. The Ti can 
now be expressed in terms of these new multipoles M.i± and £i±: 



V-^4/ 



E 

I 



( ipu, {i+m-i PI PI \ (^lA 
(/+i)p/ ipi 

_ r>ii nil 
^l+l /-I 

nil nil nil I 




PI' 





PI' 



PI' 



V £1- J 



(BIO) 



Eq. (piO|) can now be inverted to yield the electric and magnetic multipoles. Using the or- 
thogonality of the Legendre polynomials P{x) and recurrence relations to relate the deriva- 



tives P' and P" to the P's, we obtain [P,0,ia 



(2{l + l)Mi+\ 

2{l + 1)^,+ 
2l£i- 



dx 



-^dP 

21+1^-^ 

-\dP 



( Pi -Pi+i 
-Pi Pi-i 
Pi+i -Pi 

Pi wk(Pi-Pi 



1+2 

2«+3 



21+V 
{Pl+2 



\Pl 



l-l 



Pi) 

1-2) 



\ 







T2 


l+l Qp 
21+1^-^ 




dP 1 
21+1^^ I 





(Bll) 



with dP = P^+i — Pi^i. The T-matrix amplitudes T^'^ can now be obtained by multiplying 



the M.I+ and Ei± with =F Jkq /(/ + 1) (cf. Table 
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Finally, we want to give the connection of the CGLN amplitudes !Fi to the helicity 
amplitudes |1TT1| : 



Hi = -ie-<^(-i|^|l,-i) 

= -i(-i|^|l,i) 
H, = -ie-2i<A(l|^|l,-i) 

H, = -ie-'^(i|^|l,i) 

The numbers in the brackets denote the helicities of the photon and the initial/final baryon. 
In the cm. system, the spins of the nucleon are opposite to the helicities. 



--^ sin ^ cos ^{J^s + ^4) 



V2 



cos ■ 



(J^2-J^i) + sin2|(JS-J^4) 



^sin6'sin|(J^3 - Ta) 



(B12) 



V2 



cos ■ 



(^2 + J^l) +COs2|(J^^3 + J^^4 



3. Compton scattering 

In the case of Compton scattering, the decompositions of and fi analogous to 
and (P38| ) are rather lengthy and will not be given here. They can be found, for example, in 
5^. Its is more convenient to start from a set of helicity amplitudes, analogous to ([B12[) 



*2 = <^>i-i = 8^(-l,-|mi'|)' <^5 = $3 3 = 3^(l,-||T|l,-i), (B13) 

*3 = <J>|-3 = 8^(-l,|mi'i)' "^6 = $3,3 = 8^(-l,||T|l,-i). 



Again, the numbers in the brackets denote the helicities of the initial and final photon and 

fi 



baryon. Their expansion in terms of the amplitudes fif is given by 



1 ) 



$3 = iE(/ + i)v'/(/ + 2) [(/ + 2)(/(;+^)-T/S^ 



(B14) 



n = m^ + mi + 2) [(/S^)- ± f!^-) 
I 

fl+ I fl+ \lj'+2 



'^{IeE i /mm) =1= '^{fsM =1= /me)]'^3 

2 '-^ 2 

with = = Z^"]^ = /^/e '==■ 0. In an abbreviated notation, we can write ( P14|) as: 



^ss' = E(2J+ l)<^i'Cii„ with ^i, = y dx $,,,(a;)di,(x). (B15) 
J -1 

The inversion of (pi4]) can now be shown to be ||55|| : 
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fl+ 

J EE 

MM 

/■('+!)- 
J EE 



1 



MM 

1+ 

EM 

ME 

ME 
EM 



WW 
1 

/ 



1+ 

EM ) 

ME 



I 

y/l(T+2) 



4 



2(«+2) 



(B16) 



where the last hne follows from time-reversal symmetry. The T-matrix amplitudes Ti± are 
now given by: 



EE 
MM 



k/(/ + l)/ 



i± 

EE 1 
MM 



EM 
T^ME 
^1+ 



k J/(Z + 2)(/ + l) / 



/+ 

EM ■ 

ME 



(B17) 



APPENDIX C: ISOSPIN DECOMPOSITION 
1. Meson nucleon scattering 



For the 1 = 1 mesons tt and C, we use the usual projection operators with the matrix 
elements (a, b = tt, 

{hj\Pi/2\ai) = ^TjTi 

{bj\P3/2\ai) = 6ji - KjTi (CI) 

in a cartesian basis. With the help of this, all possible reactions can be written as: 

{b,\TM = \T,r,Tli' + (5,, - \r,n)Tli\ (C2) 

For the pure I = \ reactions involving vr and C the projector is usually taken to be P1/2 = t 
56| . This choice has the disadvantage that it does not agree with the Clebsch-Gordan 



coefficients for the different reactions channels. Therefore we choose instead (a = vr, 
h = Tj^k): 

{h\TM = -^nTli\ (C3) 

This has no influence on the calculated quantities, since in the end, we convert our amplitudes 
to the normal convention. 



2. Photoproduction 

In our calculation we use a decomposition slightly different from the one normally used 
in pion photoproduction (^. For the same reason as in the hadronic case, we always use 
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— -^T instead of r as projection operator for transition were either the photon or the meson 
have isospin / = 0. For the 7 = 1 mesons this leads to: 

(^'.•17^^717) = -^^K + l^^-'^T^i' + ('^.•3 - ^r,r3)T,f , (C4) 

while for the other mesons we have: 

{b\T^\l)=T^,-^T,T^/;. (C5) 

This choice guarantees that there are no additional factors to be taken into account in the 
calculation of the hadronic rescattering. 



3. Compton scattering 



For Compton scattering no isospin decomposition is performed. As has been discussed 
in Sec. PH in this case the rescattering takes place through the physical intermediate states 
(tt^p, 7r°n, etc.). Therefore, the Compton amplitudes are calculated for both channels 7p 
and ■jn directly. 



APPENDIX D: OBSERVABLES 



In this section we want to summarize the relations of the various amplitudes from App. 
to the observables. Again, more detailed accounts can be found in various places P|, ^ , |53| . 
In the following, the notation O should indicate that the amplitude O for a specific reaction 
channel has been constructed from the isospin-decomposed amplitudes: 

d = Y.p'o'. (Di) 

The factors are Clebsch-Gordan coefficients and can be determined from isospin decom- 
positions listed in App. |C[ 



1. Meson nucleon scattering 



The observables consist of the total cross sections a, the differential cross sections ^ 



and the final-state polarizations V IT 



cr 



da 

dn 



(^ + 1) 



T, 



1+ 



+ 1 



'i I 

5( = - sin 6* ^ \fi+ - fi 



Pi 



P! 



da 



\fr + \9\, ^V = -2lmirg). 



(D2) 
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2. Photoproduction 



In terms of the helicity amplitudes Hi of Eq. 



the differential cross sections and 



the three single-polarization observables can be written as [|TT 

4 



q 



— E 
dVL 

— V 

dn 

do_^ 

dn 



Re{HiHl - H2H;) 
^lm{H,H; + H2HI) 

k 



Photon asymmetry, 
Recoil polarization, 
Target asymmetry. 



(D3) 



In complete analogy to the hadronic case, the total cross sections are given in terms of the 
partial waves by: 



(!+!)( 



^1+ 



+ 



+ 



M 



(D4) 



From wich the reduced cross section can be calculated via: 



^red 



I ak 
47rq 



(D5) 



3. Compton scattering 



Again, as in the case of photproduction, the observables are best expressed in terms of 
the helicity amplitudes ( [B14|) ||55[| : 



da 

dn 

— E 
dQ 

d<7 ^ 
— V 

dn 



$1 



+ 



+ 2( 



+ 



4: 5 



- Re (($1 + $5)^3 + - $6)*: 

- Im (($1 + $5)$: - ($2 - $6)$3) = 
Finally, the total Compton cross section is given by: 



(^ + 1)( 



T, 



MM 
1+ 



+ 



EE 
1+ 



+ 2( 



ME 



+ 



T,Tl) + l{ 



nEM 



(D6) 



nME 



(D7) 
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TABLES 



Multipol 


Partial- 
welle 


J 


p 


l,p, CG 
P^ = - 


LN-Notation 

P^ 


= + 


Mi+ 

£i+ 
£i- 




i + h 




l 
I 

l + l 
l-l 


Mi+ 


l + l 
I - 1 

/ 

I 


Mii+i)- 
Ei- 



TABLE L Connection of the multipoles M.i± and £i± to the partial waves and usual CGLN 
multipoles. I denotes the photon, J the total and 1^, the meson angular momentum; a = , ^ 

^ F , ^ B , y'kq 1(1+1) 









X^JDF 


X^,/DF 


X^^/DF 


XL/DF 


X^c/DF 


X^,/DF 


X^^/DF 


SM95-pt-l 


9.61 


7.15 


13.08 


6.09 


5.17 


3.13 


5.86 


1.73 


3.28 


SM95-pt-2 


7.76 


5.20 


9.62 


3.00 


3.91 


5.78 


9.43 


1.95 


3.77 


SM95-pt-3 


5.58 


3.40 


6.69 


2.78 


4.09 


3.88 


8.10 


1.86 


4.21 



TABLE II. pci' datapoint values for the different fits. First line: Ohta's method with 
fixed hadronic parameters, second: Ohta's method with a refit of the hadronic parameters, third: 
Haberzettl's method, all parameters refitted. Also the x^/DF-values for the different reaction 
channels are given separately. 



g 


Value 


g 


Value 


K 


Value 


g 


Value 


g-KNN 


13.05 


gpNN 


2.35 


KpNN 


2.26 




1.12 




13.05 




1.94 




3.44 




1.08 




13.09 




2.07 




3.01 




1.05 


gr)NN 


1.13 


gaoNN 


0.18 






g-yrju! 


0.20 




1.01 




0.52 








0.23 




0.36 




0.41 








0.36 


gxNA 


-6.12 


gK*NA 


-6.52 


KK*NA 


-0.43 








-6.25 




-5.82 




-0.42 








-8.65 




-5.99 




-0.45 







TABLE III. Couplings of the mesons to the nucleon and g^j^pO^^ as obtained in the fits. First 
line: SM95-pt-l, second: SM95-pt-2, third: SM95-pt-3. 
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1 44^ 
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1 47Q 
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TABLE IV. Extracted resonance parameters. Notation as in Table III. the coupling gr^NR is 
given instead of the partial width, width in keV. The signs of the couplings are given in brackets. 
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TABLE V. Fitted hadronic z-parameters of the spin-| resonances. Notation as in Table 



III. 
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TABLE VI. Extracted helicity couplings (in units of 10 ^ GeV ^/^) and z-parameters for the 
l = \ resonances. The first value denotes the AP and the second A". First line: PDG- Values |17|, 
second: Arndt et al. [§, third: SM95-pt-l, fourth: SM95-pt-2, fifth: SM95-pt-3. 
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TABLE VII. Extracted helicity couplings (in units of 10 

3 



GeV ^^'^) and z-parameters for the 



1=2 resonances. Notation as in Table |V] . 
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Value 
[GeV] 




Value 
[GeV] 




Value 
[GeV] 




Value 
[GeV] 




Value 
[GeV] 


An 


1.23 
1.19 

1.24 


\h 

^h/2 


1.24 
1.24 
1.72 


^^1/2 


0.92 
0.93 
1.06 


\h 

^h/2 


1.06 
1.10 
1.13 


A^ 


3.98 
3.00 
3.59 


At 


0.70 
0.72 
0.70 



TABLE VIII. Values of the fitted cutoff-parameters A. Notation as in Table III. For the 
nucleon resonances the indices h and e name the Cutoffs at the hadronic and electromagnetic 
vertices. 
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TABLE IX. Comparison of SM95-pt-2 and a fit using PS-coupling at the i^ATA-vertex. Shown 
are only the various -values and the coupling constants gxNA, 9k* N A and kk*na- 



FIGURES 




FIG. 1. Results for the proton multipoles of pion photoproduction together with the data 

from SP97 [|. Shown are all three fits: (•••) SM95-pt-l, (- ••) SM95-pt-2, ( ) SM95-pt-3. For 

comparison we also show the results of a T-matrix calculation [^2| (— — ). 
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FIG. 4. Results for the total 7p rjp (upper plot) and jp K^A (lower) cross sections 



compared to data from pj5| , |40| . In the case of eta photoproduction the reduced cross section 
shown. The linecodes are those of Fig. |l|. 



1.491 GeV 

1.0 
0.8 



0.6^ i i i i i i i 



0.0- 



1.636 GeV 



1.711 GeV 



0.8 
0.6 
0.4 

0.2^ 



1.496 GeV 




1.644 GeV 




1 .792 GeV 



1 .502 GeV 



1.666 GeV 



1.868 GeV 




r--- """-^^ i~T-¥\ 



1.689 GeV 



1.897 GeV 



-1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 1.0 

COS 



FIG. 5. Comparison with data for the calculated differential jp 
energies. The datapoints are taken from: [04l (•), [5^, the linecodes 



ijp cross sections at different 
are those of Fig. ||. 
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FIG. 6. Polarization observables S, V and T for 7p r]p as compared to the data for different 
energies. For the photon asymmetry S also the calculation of Knochlein et al. is shown (— — ). 
The datapoints are taken from |35,38|, the linecodes are those of Fig. |l[ 
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FIG. 7. Neutron/Proton ratios and the isospin asymmetry X for 7p — > r/p as compared to the 
data for different energies. The datapoints are taken from |37], the hnecodes are those of Fig. ^ 
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FIG. 8. Comparison with data for the calculated differential 7p K^A cross sections and 

A-polarizations for different energies. The datapoints are taken from: ('—')) |40|, the linecodes 
are those of Fig. |l|. 
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FIG. 9. Comparison with data from [^^jS^ for the calculated 7p jp differential cross sections 
(top) and photon asymmetries (bottom). The linecodes are those of Fig. |l[ Also shown are the 
result of the isobar- model of Wada et al. |^l|] (— •) and the dispersion relation calculation from 
LVOV II (- -). 



1 I 40 




Vs [GeV] Vs [GeV] 

FIG. 10. Influence of the rescattering on some of the multipoles. Shown are the results for the 

fit with fixed hadronic parameters with ( ) and without (— — ) direct resonance couphng. (a): 

P33(1232) and and Mf(^-multipoles. (b): L>i3(1520) and El^_- and M2"_-multipoles. 
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FIG. 11. Comparison of contributions to the differential 7p cross section (data from |3T|). 

Full calculation using SM95-pt-2 ( ), without the vr'' and r] contributions in the t-channel (•••), 

Born s- und u-channel diagrams only (— — ). Left: ^33(1232) contributions only (—••). Right: 
1)13(1520) contribution only (—••). 
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FIG. 12. Influence of the rescattering on the nucleon s- and u-channel contributions to Compton 
scattering. Shown are the differential cross sections for two energies for the /C-matrix calculation 
( ) and using the T-matrix approximation (— — ) . 
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FIG. 13. Different contributions to Re(Eg^). Full calculation using SM95-pt-2 (- 



SM95-pt-2, but without the 5*11(1535) (— — ), (— ••) and (•••) show the results of [1^, both with 
and without the S'ii(1535). 
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FIG. 14. Sensitivity of the Compton cross section under 90° on the hehcity couphngs of the 

i:'i3(1520). Full calculation using SM95-pt-2 ( ), result using SM95-pt-2 and PDG- values for 

^1/2,3/2 °f t'^^ 1)13(1520) ( ) and the ^33(1232) contribution only for A^^^ = 4.0 GeV (- ••) and 

1.1 GeV (•••). The data from |^l[ are shown without their errorbars. The marking on the energy 
axis indicates the range of data used in the fits. 




